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We present our preliminary studies of an eflective model of a superconductor with short coher- 
ence length involving magnetic interactions. The Hamiltonian considered consists of (i) the effective 
on-site interaction U, (ii) the intersite magnetic exchange interactions (J^, J^") between nearest- 
neighbors and (iii) the intersite charge exchange term I, determining the hopping of electron pairs 
between nearest-neighbor sites. In the analysis of the phase diagrams and thermodynamic proper- 
ties of this model for half-filling (n = 1) we have adopted the variational approach, which treats the 
on-site interaction term exactly and the intersite interactions within the mean-field approximation. 
One finds that the system considered can exhibit very interesting multicritical behaviors (includ- 
ing tricritical, critical-end and bicritical points) caused by the competition between magnetism and 
superconductivity, even for n — 1. Our investigations show that, depending on the values of interac- 
tion parameters, the system at half-filling can exhibit three homogeneous phases: superconducting 
(SS), (anti-)ferromagnetic (F) and nonordered (NO). The transitions between ordered phases (SS, 
F) and the NO phase can be first order as well as second order ones, whereas SS-F transition is first 
order one. Temperature dependencies of the order parameters and thermodynamic properties of 
the system at the sequence of transitions: SS— >F— >NO with increasing temperature for J/I = 0.3, 
U /lo = 0.69 and n = 1 are also presented. 

PACS numbers: 71.10.Fd, 71.10.-w, 74.20.-z, 74.81.-g, 75.30. Fv 



I. INTRODUCTION 

There has been much interest in superconductivity 
with very short coherence length. This interest is due 
to its possible relevance to high temperature supercon- 
ductors (the cuprates, doped bismuthates, fuUerenes and 
iron-based) and also to the several other exotic super- 
conducting materials (for a review, see Refs. [T]and[l]and 
references therein). It can also give relevant insight into 
a behavior of strongly bounded fermion pairs on the op- 
tical lattices. 

The interplay and competition between superconduc- 
tivity and magnetic orderings is currently under intense 
investigations (among others in iron chalcogenides and 
cuprates, e. g. Refs. |3H5]and references therein). A con- 
ceptually simple model for studying that competition will 
be studied in this report. 

The Hamiltonian considered has the following form: 

where ni = X;<,n^a, = c+c,^, pf = {p' ^ c+cf^, 
= (l/2)(ni^ — ni^). Ci„ and denote annihilation 
and creation operators of an electron with spin a =ti i 
at the site i, which satisfy canonical anticommutation 
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relations 

{Cia,C+„-} = 5,j5„^,, {C.,„, Cjo-'} = {c+ , c+^,} = 0, (2) 

where 5ij is the Kronecker delta. indicates the sum 

over nearest-neighbor sites i and j independently. U is 
the on-site density interaction. / is the intersite charge 
exchange interaction between nearest neighbors and J is 
the Ising-like magnetic interaction between nearest neigh- 
bors, p is the chemical potential, depending on the con- 
centration of electrons: 

i 

with < rt < 2 and N is the total number of lattice sites. 

There are two competitive interaction parameters of 
the model: (i) the pair hopping interaction /, determin- 
ing the electron pair mobility and responsible for the 
long-range superconducting order (local pairing mecha- 
nism) and (ii) the Ising-like interaction J between near- 
est neighbors responsible for magnetism in the system. 
The on-site density-density interaction U contributes (to- 
gether with /) to the pair binding energy by reducing 
{U > 0) or enhancing {U < 0) its value. Moreover, re- 
pulsive U > favors magnetic ordering. To simplify 
our analysis we do not include in Hamiltonian ([T|) the 
single electron hopping term j tijcf^cja) as well as 
other inter-site interaction terms. This assumption cor- 
responds to the situation when single particle mobility 
is much smaller than the pair mobility and can be ne- 
glected. 

The interactions U, I and J will be treated as the 
effective ones and will be assumed to include all the pos- 
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sible contributions and renormalizations like those com- 
ing from the strong electron-phonon coupHng or from the 
coupling between electrons and other electronic subsys- 
tems in solid or chemical complexes [T]. 

Ferromagnetic XY-order of pseudospins pi (for / > 0) 
corresponds to the SS phase (s-pairing superconduct- 
ing), whereas the antiferromagnetic XY-order (for / < 0) 
- to the S77 phase (77-pairing superconducting). For 
tij = there is a well known isomorphism between the 
SS and Sry cases (with an obvious redefinition of the 
order parameter: A = A55 = (p"), for / > and 

A^S = j^Yji 6xp (iQ • Ri){p'i), for / < 0, Q is half of the 
smallest reciprocal lattice vector) for lattices consisting of 
two interpenetrating sublattices such as for example SC 
or BCC lattices. One should also notice that, in the ab- 
sence of the single electron hopping term, ferromagnetic 
(J > 0) interactions are simply mapped onto the antifer- 
romagnetic cases ( J < 0) by redefining the spin direction 
on one sublattice in lattices decomposed into two inter- 
penetrating sublattices. Thus, we restrict ourselves to 
the case of / > and J > 0. 

We have performed extensive study of the phase di- 
agrams of model ([TJ for arbitrary n and /i [6|. In the 
analysis we have adopted a variational approach (VA), 
which treats the onsite interaction term (U) exactly and 
the intersite interactions (/, J) within the mean-field ap- 
proximation (MFA). In this paper we present our prelim- 
inary results for the half-filling (n = 1). 

Model ([ij has been analyzed within VA only for par- 
ticular cases: (i) J = HHII] and (ii) / = [HI ES] till 
now. The rigorous results for / = in ground state have 
been also obtained [15 . Some preliminary study of the 
/ = case in finite temperatures using Monte Carlo sim- 
ulations has also been done for a square lattice pJ3, . 

The ferromagnetic (F) phase is characterized by 
nonzero value of the magnetic order parameter (magneti- 
zation) defined as m = (1/^) X^i (•^f) (and A = 0), in the 
superconducting (SS) phase the order parameter A 7^ 
(and TO = 0) and in the nonordered (NO) phase m ~ Q 
and A = 0. 

Within the VA the intersite interactions are decou- 
pled within the MFA, what let us find a grand canon- 
ical potential per site w(/i) (or free energy per site 
/(n) = i^{p) + pn) in the grand canonical ensemble. One 
can also calculate the averages: n, A and to, what gives 
a set of three non-linear self-consistent equations (for ho- 
mogeneous phases). This set for T > is solved numer- 
ically and one obtains A, to, and n (or /i) when /i (or 
n) is fixed. It is important to find a solution correspond- 
ing to the lowest a;(/x) (or f{n)). For n = 1 one obtains 
f.1 = U/2 and two equations for A and m need to be solved 
numerically. 

We also introduce the following denotation: Iq = zl, 
Jo = zj, where z is the number of nearest neighbors. 
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FIG. 1. (Color online) Phase diagrams ksT/Io vs. U / lo 
at half-fiUing (n = 1) for (a) J/I = 0.3, (b) J/I = 0.51 and 
(c) J/I = Dotted and solid lines indicate first order and 
second order boundaries, respectively. T, E and B denote 
tricritical, critical-end and bicritical points, respectively. 



II. RESULTS AND DISCUSSION 

There are two well defined limits of model ([iJ: 

(i) [/ — > — cx) favoring superconductivity and 

(ii) U -foo, where only magnetic orderings can 
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appear in the system. 

For U — >■ —oo (states with single occupancy are ex- 
cluded and only local pairs can exists in the system) 
the model is equivalent with the hard-core charged bo- 
son model on the lattice [HI HSl IS]- In this limit the 
SS-NO transition is second order one and is to the NO 
phase being a state of dynamically disordered local pairs. 
The SS-NO transition temperatures increase monotoni- 
cally with decreasing |n — 1|. The maximum value of the 
transition temperature is fesT/Io = 1 for n = 1 |9l ITT]. 

In the opposite limit (i. e. [/ — > +oo) the double oc- 
cupied sites are excluded and only the magnetic states 
can occur on the phase diagram |12l ITS] . At sufficiently 
low temperatures the homogeneous phases are not states 
with the lowest free energy and the PS state are stable (if 
n 7^ 1). On the phase diagram there is a second order line 
at high temperatures, separating the F and NO phases, 
whereas first order transition takes place at lower tem- 
peratures, leading to a phase separation of the F and NO 
phases. The critical point for the phase separation (tri- 
critical point) lies on the second order F-NO line and it 
is located at fcsT/Jo = 1/3 and n = 1/3 [IJ. The F-NO 
(second order) transition temperature decreases with in- 
creasing jn — 1| and its maximum value is fc^T/Jo = 1 
for n = 1 [13,. 



A. The phase diagrams at half-filling 

A few representative fc^T/Jo vs. U/Iq phase diagrams 
of model ([T|) evaluated for various ratios of J /I at half- 
filling (n — 1) are presented in Fig. IT] 

The phase diagram ksT /Iq vs. U /Iq for J/I = 0.3 and 
n = 1 is shown in Fig. [iji,. Two ordered phases: the 
SS phase and the F phase are separated by first order 
boundary on the diagram. Both order parameters change 
discontinuous at the SS-F transition. With increasing 
U /Iq the SS-NO transition temperature decreases from 
ksT/Io = 1 at U/Io -oo. At U/Iq = | ln(2) ~ 0.462 
and ksT/Io = 1/3 the transition changes its type from 
second order one to first order one resulting in the tri- 
critical point T on the phase diagram. The F-NO tran- 
sition temperature is slightly dependent on U /Iq and in- 
creases to ksT/Io = 0.3 {ksT/jQ = 1) at U/Iq -fcx). 
The F-NO second order line ends at critical-end point E 
on the first order boundary of the SS phase occurrence. 

The possible sequences of transitions with increasing 
temperatures and the transition orders of them are listed 
below (for J/I = 0.3): 

(i) SS-^-NO: second order, for U/Iq < 0.46 and first 
order, for 0.46 < U/Io < 0.63, 

(ii) SS— >-F— 5-NO: first order and second order, respec- 
tively, for 0.63 < U/Io < 0.7, 

(iii) F^NO: second order, for U/Ia > 0.7. 

The phase diagram for J/I ~ 0.51 is qualitatively dif- 
ferent than that for J/I — 0.3. For J/I = 0.51 the sys- 



(a) 



0.5 



0.4 



0.3 



0.2 



0.1 



0.0 









SS 


F 


NO 



0.0 0.1 0.2 0.3 0.4 0.5 

k T/I 

B 




k T/I 

B 



FIG. 2. (Color online) Temperature dependence of (a) su- 
perconducting order parameter |A| and (b) magnetic order 
parameter m for J/I = 0.3, U/Iq = 0.69 and n = 1. 



tem exhibits bicritical behavior (Fig. [T]d) in contrary to 
the tricritical behavior (and occurrence of £^-point) for 
J/I = 0.3. Similarly as for J/I = 0.3, the SS-F tran- 
sition is first order one while the SS-NO and F-NO 
transitions are second order ones. The two second or- 
der boundaries and the first order boundary merge at 
bicritical point B. 

The system exhibits the tricritical behavior for 
J/I < 0.5, whereas the bicritical behavior occurs for 
0.5 < J/I <2. For J/I > 2 the system exhibits tricriti- 
cal behavior again, however the tricritical point T is lo- 
cated at the F-NO line at U/Jq = -| ln(2) ~ -0.462 and 
ksT/Jo = 1/3 (cf. Fig.jl]:). For J/I > 2 the F-NO tran- 
sition can be second order (for U/Jq > —0.46) as well as 
first order (for U/Jq < —0.46). Notice that the axis in 
Fig. ^ are normalized by Iq, not by Jq. 

One should notice that, for any J/I, with increas- 
ing U/Iq the SS-NO transition temperature decreases 
monotonically from kgT/lQ = 1 at J7 — >■ — oo, whereas 
the F-NO transition temperature is an increasing func- 
tion of U /Iq (to its maximum ksT / Jq = 1 at ?7 — >■ -|-cx)). 

Let us concentrate now on temperature dependencies 
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of the order parameters and thermodynamic properties 
of the system at the sequence of transitions: SS^F— )-N0 
for J/ 1 = 0.3, U/Iq = 0.69 and n = I. 



B. The order parameters 

The temperature dependencies of the order parame- 
ters: A and m for J /I ~ 0.3, U /Iq — 0.69 and n ^ 1 are 
presented in Fig. [2j It is clearly seen that at the SS-F 
transition (at kBTd/Io = 0.16) the both order parame- 
ters: superconducting order parameter A and magneti- 
zation TO change discontinuously. In the SS phase A ^ 
and TO = whereas in the F phase to ^ and A = 0. The 
F-NO transition (at kBTc2/Io = 0.24) is connected with 
a continuous decay of to at the transition temperature. 



C. The thermodynamic properties 

Calculating the free energy per site / one can obtain 
thermodynamic characteristics of the system for arbi- 
trary temperature. The double occupancy per site D 
is defined as: 



D 



N ^ 



dU 



(4) 



and it is related with the local magnetic moment 7 by 
the following formula: 



i i 



(5) 



df 



-n- D, 



because \h^^ - ^4! = ('^it - "-4)^ = ^il + fiii - 2%f^4, 
nf^ = fiia = 0, 1 and In^^ - ri^l = 0, 1. 

The entropy s and the specific heat c can be derived 
as: 



c=-T 



(6) 



The temperature dependencies of the thermodynamic 
parameters for J// = 0.3, C///o = 0.69 and n^l are 
shown in Fig. [3] 

The concentration of paired electrons Up = 2D (nor- 
malized to the total electron concentration n) as a func- 
tion of temperature is presented in Fig. |3^. At the SS-F 
transition large amount of electron pairs is destroyed. 
Thus Up has a sharp break at the SS-F transition tem- 
perature Tci and a substantial fraction of single particles 
exists above Td- As temperature is lowered, the conden- 
sate growths both from a condensation of pre-existing 
pairs and from binding and condensation of single parti- 
cles. At the F-NO transition (at Tc2) rip is continuous. 
In the NO phase it increases to rip —)■ 0.5 at T — )■ -l-oo 
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FIG. 3. (Color online) Thermodynamic parameters (a) the 
concentration of paired electrons Up/n = 2D/n, (b) the en- 
tropy s and (c) the specific heat c as a function of fesT/Zo for 
J/I = 0.3, U/Iq = 0.69 and n = 1. 



(two electrons at the site is one of four equal probable 
configurations at the site and n = (hi) = 1). 

The temperature dependencies of the entropy s and 
the specific heat c are shown in Figs. [SJd and|3]3, respec- 
tively, s increases monotonically with increasing T. At 
Tci the entropy s is discontinuous whereas it is continu- 
ous at Tc2- One can notice that in the high-temperature 
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limit the entropy s/ks ^ lii(4) w 1.386 (there are four 
possible configurations at each site). The peak in c{T) is 
associated with the first order transition (at T^), while 
the A-point behavior is typical for the second order tran- 
sition (at Tc2). 



III. FINAL REMARKS 

We have studied a simple model of a magnetic super- 
conductor with very short coherence length (i. e. with the 
pair size being of the order of the radius of an effective 
lattice site) and considered the situation where the single 
particle mobility is much smaller than the pair mobility 
and can be neglected. 

One has found that the system considered for n = 1 ex- 
hibits various multicritical behaviors (determined by the 
ratio J / 1) including tricritical, critical-end and bicritical 
points. It has been shown that, depending on the values 
of interaction parameters, three homogeneous phases: su- 
perconducting, (anti-)ferromagnetic and nonordered oc- 
cur on the phase diagrams of model ^ at half-filling. 
The transitions between ordered phases (SS, F) and the 
NO phase can be first order as well as second order ones, 
whereas the SS-F transition is first order one. For n ^ 1 
several types of phase separated states could be also sta- 
ble in definite ranges of model parameters [6j . 

The other result of the interplay between magnetism 
and superconductivity could be appearance of triplet 
pairing [TT]. Such a solution could appear together 
with ferromagnetic spin ordering. However, in model ([T]) 
which assumes tij = such a state cannot be found. To 
investigate the possibility of occurrence of a supercon- 
ducting state with triplet pairing, the model should be 



extended to the case of finite bandwidth {tij ^ 0) and 
be analyzed taking into account intersite pairing (in par- 
ticular triplet pairing), e. g. using Hartree-Fock broken 
symmetry framework [rT-TFI . 

The mean-field approximation used to the intersite 
term is best justified if the lij and Jij interactions are 
long-ranged or if the number of nearest neighbors is rela- 
tively large. The derived VA results are exact in the limit 
of infinite dimensions d — > +oo, where the MFA treat- 
ment of the intersite interactions / and J terms becomes 
the rigorous one. 

Let us point out that in the MFA, which does not 
take into account collective excitations, one obtains the 
same results for the U-I-J^ model, i. e. model (jlJ, 
and the U-I-J^^ model, where the term s| is re- 

placed with i^t^J + ^^^7)' describing interactions 
between xy-components of spins at neighboring sites, 
sf = cf^Cii = {s~y ■ In both cases the self-consistent 
equations have the same form, only a magnetization 
along the z-axis becomes a magnetization in the xy- 
plane [12 . 
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